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Abstract 

An experimental study was conducted to investigate the structural stability of ethylene glycol-based titaniunn 
dioxide nanoparticle suspensions (nanofluids) prepared by two-step method. The effects of particle concentration, 
fluid temperature, shear rate and shear duration were examined. Particle size and thermal conductivity 
measurements in quiescent state indicated the existence of aggregates and that they were stable in temperatures 
up to 60°C. Shear stability tests suggested that the structure of nanoparticle aggregates was stable in a shear 
interval of 500-3000 s'^ measured over a temperature range of 20-60°C. These findings show directions to resolve 
controversies surrounding the underlying mechanisms of thermal conduction and convective heat transfer of 
nanofluids. 



Introduction 

Nanofluids are suspensions of nano-sized particles in 
liquids, where particle sizes are preferably below 100 
nm. At modest particle concentrations, the thermal con- 
ductivity, forced convective heat transfer, and critical 
heat flux of nanofluids were reported to be superior to 
respective base liquids [1-8]. In the backdrop of conven- 
tional heat transfer technologies approaching their 
upper limits, nanofluids are seen as a potential conten- 
der for small- and large-scale thermal applications 
[9-12]. A number of attempts had been made in the 
past, and postulates were put forward to explain the 
underlying mechanisms. Although yet inconclusive, the 
nanoparticle aggregation in liquids is believed to be one 
of the principal mechanisms behind the enhanced ther- 
mal conductivity and convective heat transfer [13-16]. 
In either case, the importance of particle aggregation 
and their stability were underlined. 

On the other hand, the aggregation of nanoparticles is 
found to be the key mechanism behind the increase of 
nanofluid viscosity and shear thinning behaviour 
[14,17,18]. Recently, it was shown that the high shear 
viscosity of nanofluids could accurately be predicted by 
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combining the conventional Krieger and Dougherty 
model and aggregation effects [18-20]. Those postulates 
were based on the assumption that, in the shear flow 
field, the aggregates will be stable because the hydrody- 
namic forces are insufficient to break the aggregates 
down to primary particles. However, the experimental 
evidences are insufficient to showcase the stability and 
particle structuring of nanofluids in flow conditions. 

In the present study, the ethylene glycol (EG)-based 
Titania (Ti02) suspensions are selected to investigate 
the stability of nanofluids in quiescent and shear flow 
fields. Also their thermal conductivities are measured at 
various temperatures and compared with theoretical 
predictions. The experimental conditions were chosen 
resembling the possible industrial applications for nano- 
fluids. Considering the bounded yet deep focus of the 
stability of nanofluids under different conditions, this 
article is reported as a letter without comparing the data 
with the other literature. 

Experimental 

Nanofluids were formulated using Ti02 nanopowder and 
EG. The dry Ti02 nanopowder purchased from Degussa 
Corporation in Germany was claimed to be consisting of 
spherical particles of 25-nm diameter. Electron micro- 
scopy (EM) imaging such as in Figure 1 suggests that the 
particles were in the form of agglomerates. In order to 
manufacture a stable nanoparticle suspension, a sequence 
of processes were followed. Further details of formulation 
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can be found elsewhere [7,21]. The EM images of the 
nanofluid confirmed that the nanoparticles were well dis- 
persed. Moreover, the light-scattering data collected 
using the Malvern Zetasizer-nano device showed that the 
suspended particles were in the order of around 130 nm 
in size. This is an indication of the formulation technique 
substantially reducing the aggregate size but failing to 
break them down to primary particles. This observation 
agrees with the recently concluded International Nano- 
fluids Property Benchmarking Exercise (INPBE) [22]. 
These nanofluids were stable for 2 months without a visi- 
ble separation, indicating the stability of aggregates in the 
long run. 

Measurements of thermal conductivity (/c, W/mK) of 
Ti02-EG nanofluids were conducted using the state-of- 
the art Lambda meter device acquired from PSL Mea- 
surement Systems GMBH of Germany. This instrument 
works on transient hot wire principle. For calibration 
with EG, the instrument reproduced the data up to 99% 
precision. 

Shear flow field was applied to the samples using a 
Bohlin rotational rheometer. The experimental condi- 
tions were as follows: shear rates 500, 1000, 2000 and 
3000 s'^; time durations 5, 10, 20 and 40 min; and tem- 
peratures 20, 30, 40, 50 and 60°C. These temperature 
and flow parameters were so chosen to suit possible 
industrial applications [19]. The shearing was preceded 
and followed by particle size measurements using Mal- 
vern Zetasizer-nano. The size measurements were 
repeated six times, and the reproducibility of data fell 
within error of 4%. In all instruments, the thermal equi- 
librium was ensured by leaving the samples at measur- 
ing temperature for a sufficient period of time before 
taking the readings. 



Results and discussion 

Thermal conductivity (/c, W/mK) data for the samples 
are presented in Figure 2. The trends of k of the nano- 
fluid and base liquid appear alike. This follows that the 
presence of nanoparticles at these concentrations has 
not altered the dynamics of the base liquid. Interest- 
ingly, this was the case even at 60°C, indicating quies- 
cent flow fields. Also shown in Figure 2 are the 
percentage (%) enhancements of thermal conductivity. 
At any given temperature, the enhancement has system- 
atically increased with loading. However for a given con- 
centration, the enhancement appears to be fairly stable 
with temperature. This is a trend that agrees with the 
more recent literature on this area [23,24]. Also noted 
from Figure 2 are the low particle loadings unable to 
cause noticeable enhancement. This observation contra- 
dicts a section of the old literature, while agreeing with 
majority of recent study including INPBE [22] partici- 
pated by dozens of nanofluids research institutions. 

The average enhancement for each concentration in 
Figure 2b is plotted in Figure 3 together with the predic- 
tions of classical Hamilton-Grosser (H-C) model based 
on well-dispersed particles [25] and modified H-C 
model [20] based on aggregated particles. The classical 
H-C model can be written as 



fep + (n - l)ko - (n - l)(p[ko - fep) 
fep + (n - l)ko + (p{ko - fep) 



(1) 



where /c, ko, kp are, respectively, the thermal conduc- 
tivities of the nanofluid, base liquid, and particle mate- 
rial, and n is the shape factor given by n = 3/y/ with y/ 
the surface area-based sphericity {y/ = 1.0 for spheres). 

Modified H-C model based on aggregated particles 
takes the form of [20] 



k/ko 



h + jn- l)ko - (n - l)(p^{ko - fea) 

fea + (n - l)ko + (p^{ko - fea) 



(2) 



where /Ca is the thermal conductivity of aggregates 
which is estimated by the Bruggeman model for spheri- 
cal particles [26]: 



ko 



(3(^in-l)-;^ + (3(l-<^in)-l) + 



ko 



^(3^in-l)^ + (3(l 



(Pin) - 1) 



(3) 



Here, 0a is the effective particle volume fraction 
given by (p^i = (piaja)^'^ according with the fractal 
theory, and (p^^ is the solid volume fraction of aggre- 
gates given by (pin = {aja)^'^. Also a and are the 
radii of primary nanoparticles and aggregates, respec- 
tively [27], and D is the fractal index having a typical 
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Figure 2 Thermal conductivity of Ti02-EG nanofluids. 



value of 1.8 for nanofluids [20]. From Figure 3, the 
conventional H-C model underpredicts the measure- 
ments by a considerable margin can be seen. However, 
the modified H-C model that takes into account the 
aggregates of nanoparticles agreed well with the 
experimental data. 



Overall view of Figures 2 and 3 suggests that (i) the 
aggregation of nanoparticles is a principal mechanism 
that drives the thermal conductivity enhancement and 
(ii) the aggregates are stable in quiescent flow fields 
even at temperature as high as 60°C. Independence of 
the experimental data on temperature further suggests 
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the weak or negligible effect of particle Brownian 
motion on reported enhancement. 

Featured in Figures 4 and 5 are the studies on particle 
size in shear flow fields. All samples have the measured 
particle sizes considerably larger than the primary size 
(25 nm) reconfirming the existence of the aggregates. 
Yet, the average particle diameter {d) exhibits a narrow 
fluctuation between 126 and 132 nm, which falls within 
the boundaries of experimental error. Moreover, the 
shear rates and shear durations shown on Figure 4 had 
been unable to break the aggregates. The aggregates 
were therefore sufficiently stable under these conditions. 

Figure 5 illustrates the dependence of the measured 
particle sizes on the measuring temperature and particle 
concentration. At any given concentration, a tempera- 
ture increase of threefold (from 20 to 60°C) has not 
registered a notable size change. Here, the indication is 
the temperature stability of aggregates. Furthermore, a 
concentration increase by 16 folds (from 0.5 to 8 wt%) 
has caused only a modest increase in size which again 
falls within the experimental error. 

Conclusions 

Experiments were conducted to study the dependence of 
shear stability of nanofluids on temperature, particle 
loading and shear rate. Observed weak dependence of 
thermal conductivity enhancement on temperature sup- 
ports the claim of particle aggregation as a principal 
mechanism behind the enhancement. Moreover, the 
aggregates in quiescent flow fields were stable in tem- 
peratures up to 60°C. The data on shear stability show 
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Figure 5 Average particle size after 40 min of shearing at 3000 s \ 



that the aggregates are sufficiently stable over a range of 
rigorous shear rates and temperatures. The observations 
of thermal conductivity and particle size complement 
each other in terms of predicting the former from the 
latter. A comparison of the present findings with the lit- 
erature data is currently underway and will be reported 
in future. 
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